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We demonstrate that chiral symmetry restoration in quenched finite temperature 

O 

QCD depends crucially on the Z3 phase of the Polyakov loop V . This dependence is 



a general consequence of the coupling of the chiral order parameter to the Polyakov 
loop. We construct a model for chiral symmetry breaking and restoration which 
includes the effect of a nontrivial Polyakov loop by calculating the effective potential 

o 

fSj ' for the chiral condensate of a Nambu-Jona-Lasinio model in a uniform temperature 

dependent Aq gauge field background. Above the deconfinement temperature there 

OS ■ 



are three possible phases corresponding to the Z3 symmetric phases of the Polyakov 



loop in the pure gauge theory. In the phase in which tr c ("P) is real and positive the 
first order deconfining transition induces chiral symmetry restoration in agreement 
with simulation results. In the two phases where Re\tr c (V)] < the sign of the 
leading finite temperature correction to the effective potential is reversed from the 
normal phase, and chiral symmetry is not restored at the deconfinement transition; 
this agrees with the recent simulation studies of Chandrasekharan and Christ. In 
the case of SU(N) a rich set of possibilites emerges. The generality of the mecha- 
nism makes it likely to occur in full QCD as well; this will increase the lifetimes of 
met ast able Z3 phases. 
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I. INTRODUCTION 



Restoration of chiral symmetry is observed at high temperatures in simulations of both 
quenched and full QCD. In interesting recent work Chandrasekharan and Christ have shown 
through simulation that the chiral properties of quenched finite temperature QCD depend 
strongly on the Z 3 phase of the Polyakov loop V |IJ. Only in the phase where Re[tr c (V)] > 
and Im[tr c (V)] = is chiral symmetry restored at the deconfinement transition. In the 
other two Z3 phases, chiral symmetry is not restored. In this letter we provide a theoretical 
explanation for these results and make predictions for the general case of SU(N). 

The physical mechanism which couples the chiral condensate to the Polyakov loop is 
simple: finite temperature corrections to the quark determinant involve quark trajectories 
which are topologically non-trivial in the Euclidean temporal direction. Such trajectories 
are naturally weighted by powers of the Polyakov loop. When Re[tr c (V)} changes sign, 
the sign of the leading finite temperature correction to the effective potential for the chiral 
condensate changes as well. 

In our efforts to understand the results of Chandrasekharan and Christ, we have studied 
an effective U(Nf) x U(N C ) chirally invariant Nambu-Jona-Lasinio (NJL) model coupled 
to a uniform temperature dependent A gauge field. NJL models have been analyzed in 
great detail, including their finite temperature properties 0,0]. We extend these models by 
allowing a dependence on the uniform external gauge field Aq. The Aq influence on the phase 
structure of this model manifests itself solely through the Polyakov loop. The Polyakov loop 
is taken to be identically zero for temperatures below the deconfinement temperature Td 
and jumps to a nonzero value at Td, as in pure 577(3) gauge theory. The deconfining phase 
transition is known to be first order for SU(3) and believed to be first order for all SU(N) with 
N > 3; SU(2) exhibits a second order phase transition and must be considered separately. 

We use the effective potential formalism to study our extended model. Since closed quark 
loop effects are by definition absent in quenched QCD, it may at first glance seem strange 
to consider an effective potential which contains a contribution from the quark determinant. 
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However, consider diagram (a) in Fig. |]. It is one of an entire class of diagrams which 
presumably contribute to chiral symmetry breaking in the quenched approximation. It 
has long been argued that such diagrams are treated approximately in NJL models by the 
interaction shown in diagram (b) [[|. If we regard the effective potential as simply a tool 
to derive the more fundamental gap equation for the chiral condensate, then the potential 
must include terms such as the one shown in diagram (d) of Fig. [I], since they are obtained 
from Fierz transformations of quenched QCD terms like the one shown in diagram (c). 



II. EFFECTIVE POTENTIAL 

We begin with a Euclidean action of the form 

Nj-l 

£ E = y(p + m)V -G [(^A a ^) 2 + (^ 75 A a ^) 2 ] (2.1) 

a=0 

where the quark field \1/ carries flavor and color ||. The A a are the generators of U(Nf). 
The covariant derivative reflects the presence of a uniform background temporal SU(N C ) 
gauge field: 

D„ = + igA^(x) with A„(x) = 5 0fl A . (2.2) 

Let 5*o denote the propagator for a free quark with mass m, and let 5* denote the quark 
propagator of the interacting theory. Following Cornwall, Jackiw, and Tomboulis f|, an 
effective action r(S') for the theory in Eq. ( |2.1|) is given by 



T(S) = Tr [in (S^S)] - Tr (s^S - l) - G ^(diagrams), (2.3) 

where the diagrams are shown in Fig. [| with each internal fermion line being associated with 
the propagator S. The two rightmost diagrams sum to zero if S(x, x) is a scalar in spinor 
space. In the presence of a uniform background A gauge field this scalar property does not 
hold at finite temperature. For simplicity, however, we consider only the first diagram. This 
is the Hartree approximation and the leading term in a 1/N C expansion. It follows that 
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T(S) = Tr 



Tr (S^S -l)-2G dt d 3 x Kd [S f (x,x)]} 



(2.4) 



where the sum is over all flavors and the subscripted trace is over Dirac and color indices 
alone. (3 is the inverse temperature. The equation for the propagator is 



S f l (x,y) = S f(x,y) + AGtx cd [S f (x,x))6(x,y). 



(2.5) 



Letting 07 = tr c a [Sf(x,x)], the quark mass of a single flavor in the interacting theory is 
given by Mf = rrif + 4G<jf. Conventionally, rrif and Mf are considered to be the current 
and constituent quark masses, respectively. 

The phase structure of the model is conveniently studied using the effective potential 
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d 3 k 



tr r In 



(k + gA ) 2 + u 2 % 



(2.6) 



where u ^ = [k 2 + (mj + AGaf) 2 ] 1 ^ 2 , and the sum over k is the mode sum over Matsubara 
frequencies. Evaluating this mode sum and discarding an irrelevant constant, 

d 3 k 



V(a,A ,T) = J2 
f 



2Ga) - 2N C 



(2n 



c\ OO 

--T 



-l) n 



n=l 



n 



(2.7) 



where V = exp(i/3gA ) is the Polyakov loop associated with the background field A . The 
second term in Eq. (|2.7|) is an integral representation of the zero-temperature quark determi- 
nant which sums all possible zero-point energies. The infinite series of integrals in the third 
term represents the finite temperature corrections. Intuitively, each order in n is associated 
with quark paths that wind around space-time in the temporal direction a net number of 
times n. This is essentially an image expansion of the finite temperature quark determinant 
0J^]. The standard NJL model is recovered by setting A to zero. 

The evaluation of the quark determinant at finite temperature is slightly subtle. Af- 
ter regularization, the zero-temperature portion of the determinant contains an M 4 log(M) 
term. This logarithmic term is cancelled by a similar term from the finite temperature 
contribution; every order in the summation over n contributes to this cancellation. On the 



other hand, the integrals associated with each order in n can be performed analytically, 
with no need for the imposition of a cutoff, resulting in modified Bessel functions K n . It is 
tempting to approximate the finite temperature correction by simply keeping the first few 
terms in the series. However, this approximation fails to reproduce the correct behavior of 
the effective potential for small constituent masses and can lead to unphysical behavior. 

The image expansion nevertheless remains an effective tool for understanding the role 
of Zn phases in critical behavior. The basic physics for the change in critical behavior 
associated with the Polyakov loop is given by the n = 1 term in Eq. (|2.7|) . Terms higher 
order in n are suppressed by powers of exp(— M/T), making this a reliable approximation 
when T « M. Keeping only the n = 1 term, the effective potential can be written as 

V(a, A , T) = ^Gaj - 2N C J 

K2 (m±^y^ v+V ^. (2 . 8) 

The crucial point is that the sign of the n = 1 finite temperature correction to V(a, A ,T) 
changes with the sign of tr c (V + T>'). If the effect of finite temperature corrections is to raise 
the normal minimum relative to M = 0, phases with Re[ti c (V)] < will instead lower this 
minimum. This suppresses chiral symmetry restoration in these phases. 

A quantitative demonstation is best performed with the full effective potential. It is 
convenient to introduce a further approximation to powers of the Polyakov loop: 

tr c (n « N c . (2.9) 

This approximation neglects the formation of triality zero states in the gluon plasma. That 
is, each quark moves independently in its own fixed background. We also regulate V(a, Ao, T) 
by introducing a non-covariant cut-off, A. Denoting the trace of the Polyakov loop by 0, 
Eq. ( p.7|) now can be resummed as 

V(a,<p,T)=Y,\2GcT 2 f -^ dkk 2 u fk 

j ^ 7T J 
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+ 



(m f + 4Ga f )T 



7T 



+ 



dk k In 



7T 2 /3 Jo 



1 + e^ fc 



+ 



+ e 



(2.10) 



N c J ~ \N* t 

In this approximation, if is zero, all finite temperature effects are absent from the potential. 
Other approximations are also defensible. Consider 



+ 



7l 2 P Jo 



dk k 2 In 



1 _|_ e -/3w/* 



+ 



+ e 



(2.11) 



which differs from the previous expression only in the last term, where the Polyakov loop 
for the quark is cancelled against its inverse for the antiquark. In this approximation, some 
finite temperature effects of quark-antiquark pairs are taken into account. In fact, even 
when = 0, there is a chiral symmetry restoration transition at exactly twice the critical 
temperature of the ordinary NJL model. 

III. DISCUSSION AND NUMERICAL RESULTS 



Following Ref. ||, we determine the noncovariant cut-off A and the constituent quark 
mass M by fixing the chiral condensate and the pion decay constant at zero temperature: 
a = (246.7 MeV) 3 and f 2 = (93 MeV) 2 . For the remainder of this letter we will use the 
numerical solution M = 335 MeV, A = 631 MeV, and GA 2 = 2.2, which assumes a current 
quark mass m = 4 MeV and N f = 2 PJ^-|10[ . 



The standard NJL model is recovered by setting <p = N c in Eq. ( 2.10 ). Chiral symmetry 
is restored by a second-order phase transition; with the A and G given above, the critical 
temperature is T c = 194.6 MeV. Our model will behave as the zero-temperature solution of 
the standard NJL model until the Polyakov loop develops an expectation value; in particular, 
chiral symmetry will remain broken while = 0. Thus, in our model the critical temperature 
for chiral symmetry restoration T c will always be greater than or equal to the temperature 
for deconfinement T^. We have shown recently in Ref. [|TTJ that the unquenched version of 
this model has T c strictly greater than Td for certain ranges of phenomenological parameters. 
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However, in lattice simulations it appears that T c = Requiring this behavior places a 
mild constraint on this model as shown in Fig. [5| At this time, the constraint is not a strong 
test of the model, because is not known: the Polyakov loop is multiplicatively renormalized 
in a way that depends on the regularization scheme. In particular, lattice simulations do 
not directly measure 0. 

If the Polyakov loop is in one of the Z3 phases of QCD in which + 0* < 0, the potential 
V(<jf,<f),T) will not have a minimum at 07 = 0, unless the temperature is significantly 
greater than the T c of the corresponding real phase. This is precisely the behavior observed 
in quenched QCD simulations Q. In Fig. [|we compare the effective potentials of the real 
and complex Z 3 phases for = 2.0 and T = 252.1 MeV to the zero-temperature effective 
potential where = 0. The relative depth of the well in the complex phase has increased as 
expected given the change in sign of + 0* . 

The general case of a theory with SU(N) gauge fields can also be considered. The 
mechanism will be the same whether the full theory or a phenomenological, low energy 
approximation is studied. The sign of the most important finite temperature correction to 
the effective potential depends on the sign of the real part of 0. In general, the temperature 
dependence of the chiral condensate will be different in different phases. Only phases in 
which the Polyakov loops are related by complex conjugation will have the same behavior. 
Thus, for both SU(2N) and SU{2N + 1), there will be N + 1 distinct behaviors. 

The finite temperature behavior when 0^ is purely imaginary is determined by higher 
order terms in the image expansion [see Eq. ( |2.7| )1- In SU(4:), scaling the temperature by 
a factor of two yields a potential identical to the original NJL model, except the remaining 
exponential in Eq. ( [2.10| ) is multiplied by a factor of |0d| 2 /A^. Naively, this would imply 
that the critical temperature is more than double that of the phase where is real and 
positive. 

In the high temperature phase of full QCD, the phases are not equivalent; quark 
effects lift the Z^ degeneracy and only the phase in which the Polyakov loop is real and 
positive is stable. The other phases are metastable, at least in the standard Euclidean 
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formalism. The action of the bounce solution has been calculated under the assumption 



that chiral symmetry has been restored in all phases ||12|| . However, it seems likely that the 
same mechanism discussed here for the quenched approximation is also operative in the full 
theory. Because the strength of the symmetry breaking decreases with increasing mass, 
the absence of chiral symmetry breaking in those phases will tend to decrease the action of 
the bounce solution, and therefore increase the lifetime of the metastable states. This will 
also modify the thermodynamic properties of these metastable phases. 
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FIGURES 

FIG. 1. Feynman diagrams contributing to the constituent quark mass in quenched QCD 
(a) and the NJL model (b) and related diagrams contributing to the effective potential (c and d, 
respectively). 

FIG. 2. Feynman diagrams contributing to the effective action. 

FIG. 3. Boundary in the {4>di T^-^Aane between parameter values which yield one phase 
transition and parameter values which yield two. Only points above the curve are compatible with 
simulation results. 

FIG. 4. Comparison of the effective potentials of the real and complex Z3 phases for <p = 2.0 
and T = T x {4>) = 252.1 MeV to the zero-temperature effective potential where = 0. 
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